Phase transitions in the boson-fermion resonance model in one dimension 
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S3 ' 

' We study ID fermions with photoassociation or with a narrow Fano-Feshbach resonance described 

by the Boson-Fermion resonance model. Using the bosonization technique, we derive a low-energy 
■ Hamiltonian of the system. We show that at low energy, the order parameters for the Bose Conden- 

sation and fermion superfluidity become identical, while a spin gap and a gap against the formation 
of phase slips are formed. As a result of these gaps, charge density wave correlations decay expo- 
nentially in contrast with the phases where only bosons or only fermions are present. We find a 
Luther-Emery point where the phase slips and the spin excitations can be described in terms of 
pseudofermions. This allows us to provide closed form expressions of the density- density correla- 
i-G , tions and the spectral functions. The spectral functions of the fermions are gapped, whereas the 

(y>j ■ spectral functions of the bosons remain gapless. The application of a magnetic field results in a loss 

D ' of coherence between the bosons and the fermion and the disappearance of the gap. Changing the 

detuning has no effect on the gap until either the fermion or the boson density is reduced to zero. 
Finally, we discuss the formation of a Mott insulating state in a periodic potential. The relevance 
of our results for experiments with ultracold atomic gases subject to one- dimensional confinement 
is also discussed. 



I. INTRODUCTION 
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Since the discovery of Bose-Einstein Condensation (BEC) of atoms in optical traps, the field of ultracold atoms has 
experienced tremendous developments in the recent years. |l| A first important step has been the use of Fano-Feshbach 
resonances 0, to tune the strength of atom-atom interaction. 0, Q Fano-Feshbach resonances take place when the 
energy difference between the molecular state in the closed channel and the threshold of the two-atom continuum in 
0^ ' the open channel, known as the detuning v, is zerofff . Near a Fano-Feshbach resonance, the atom-atom scattering 
length possesses a singularity. For v > 0, atoms are stable, but the existence of the virtual molecular state results in an 
effective attraction. For v < 0, the molecules are formed and possess a weakly repulsive interaction. Since the value of 
v can be controlled by an applied magnetic field, this allows to tunc the sign and strength of the atomic and molecular 
interactions. H, I p. H I, Q(| In particular, the use of Fano-Feshbach resonances has allowed the observation of pairs of 
fermionic|lll lUlla. IT^j or bosonic fTa. ll6L UtI Il8j atoms binding together to form bosonic molecules. At sufficiently 
low temperature, for v < 0, these molecules can form a Bose-Einstein condensate. In the case of a fermionic system, 
for v > 0, due to attractive interactions a BCS superfluid is expected. Since the BEC and the BCS state break the 
same U(l) symmetry, a smooth crossover between the two states is expected as v is tuned throu gh the resonance. 
Indeed, the BEC of molecules [Hill HJ and the crossover to a strongly degenerate Fermi gas [2lH22ll23l.l24j have been 
observed as a gas of cold fermionic atoms is swept through the Fano-Feshbach resonance. Measurement of the radio- 
frequency excitation spectraj2j| and of the specific heat |2g as well as observation of vortices in a rotating system^ 
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q \ on the v > side revealed the presence of a superfluid BCS gap, thus proving the existence of a BEC-BCS crossover. 
Such a crossover is naturally described by the boson-fermion mod el. J28l l30t l3ll 13^ . 1331 133. |3j| first introduced 
. £^ \ in the 1950s in the context of the theory of suDerconductivit yl36l l37l an d later reinvestigated in the 1980s in the 
context of polaronic|38| and high-Tc superconductivity theory. [391 1401 14]| A second important parallel development 
has been the possibility to form quasi-lD condensates using anisotropic traps |42l 1431 l44t l45j . two-dimensional optical 
lattices MM, H ElM E3 or atoms on chips.^ In one dimensional systems interactions are known to lead to a rich 
physics. |53j In particular, strongly correlated states of fermions, where individual particles are replaced by collective 
spin or density excitations, are theoretically expected. |53l l54l l55| When the interactions between the fermions are 
repulsive, both the spi n and density fluctuations are gapless with linear dispersion and this state is known as the 
Luttinger liquid[53l |5(| E3- For attractive interactions between the fermions, the spin degrees of freedom develop 
a gap, yielding a state known as the Luther-Emery liquid. [EH. liSSj Similarly, bosons are expected to be found in 
a Luttinger liquid state, with individual particles being replaced by collective density excitations [53l IB"!. l59l l60| . 
Moreover, strong repulsion can lead to the fermionization of interacting bosons i.e. the density matrix becomes 
identical to that of a non-interacting spinless fermion system, the so-called Tonks- Girardeau (TG) regime. [6ll Ir3^| 
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Experiments in elongated traps have provided evidence for one-dimensional fluctuations [42l I43L |44| . However, in 
these systems, the bosons remain weakly interacting. With two-dimensional optical lattices, it is possible to explore 
a regime with stronger repulsion. In particular, it was possible to observe the TG regime with 87 Rb atoms|4g by 
increasing the transverse confinement. The TG regime can also be reached by applying a ID periodic potential along 
the tubes to increase the effective mass of the bosons 49]. Using a stronger ID potential, it is possible to drive a 
one-dimensional Mott transition between the superfluid state and an insulating state[H3|- Another characteristic of 
atoms in a one-dimensional trap is that transverse confinement can give rise to a type of Fano-Feshbach resonance as a 
function of the trapping frequency called the confinement induced resonance (CIR) . j63L l64l l65j Recently, experiments 
have been performed on 40 K fermionic atoms in a one dimensional trap forming bound states either as a result 
of Fano-Feshbach resonances or of CIR.l6a| Both types of bound states have been observed and the results can be 
described using the Boson- Fermion model. [67j This prompts the question of whether a one dimensional analogue of the 
BEC-BCS crossover could be observed in such a system. It is well known that in one dimension, no long range BEC or 
BCS order can exist. [g^lg^ However, quasi-long range superfluid order is still possible. For fermions with attractive 
interactions, it was shown using the exactly solved Gaudin-Yang modeljlilEI that for weakly attractive interactions, 
a Luther-Emery state with gapless density excitations and gapful spin excitations was formed, whereas for strongly 
attractive interactions the system would crossover to a Luttinger liquid of bosons. |72l l73| The boson- fermion model 
was also considered in the case of a broad Fano-Feshbach resonance. (7^ In that case only bosons or fermions are 
present (depending on which side of the resonance the system is) and the results are analogous to those obtained with 
the Gaudin-Yang model. In fact, in the three dimensional case, it is possible to derive a mapping of the boson-fermion 
model with a broad resonance to a model with only fermions and a two-body interaction. |75j In the narrow resonance 
case, such a mapping is valid only very close to the resonance. It was therefore interesting to investigate what happens 
in one dimension in the case of a narrow resonance. Indeed, in the latter case, it has been shown previously |76l l77j] 
that a richer phase diagram could emerge with a phase coherence between a fluid of atoms and a fluid of molecules 
at weak repulsion and a decoupling transition for stronger repulsion. Analogous effects have been discussed in the 
context of bosonic atoms with a Fano-Feshbach resonance in[2a|. Due to the concrete possibility of forming ID Fermi 
and Bose gas with optical lattices 0,H3 some of the theoretical predictions in the narrow resonance case may become 
testable experimentally in the future. Experimental signature of the phase coherence between the two fluids include 
density response and momentum distribution function. In the present paper, we investigate in more details the phase 
in which the atomic and the molecular fluid coexist. In particular, we study the equilibrium between the atomic 
and the molecular fluid as the detuning is varied. Also, we investigate the effect of placing the system in a periodic 
potential and show that the phase coherence between the atomic and molecular fluid hinders the formation of the 
Mott state in systems at commensurate filling. Such conclusion is in agreement with a study in higher dimension[79]| . 

The plan of the paper is the following. In Sec[H]we introduce the boson-fermion Hamiltonian both in the lattice 
representation and in the continuum. We discuss its thermodynamics in the limit of an infinitesimal boson-fermion 
conversion term and show under which conditions atoms and molecules can coexist. In Sec lIIII we derive the bosonized 
expression for the boson-fermion Hamiltonian valid in the region where atoms and molecules coexist. This Hamiltonian 
is valid for a system in an optical lattice provided it is at an incommensurate filling (i.e. with a number of atoms per 
site which is not integer). We show that for not too strong repulsion in the system, a phase where the atomic and the 
molecular superfluid become coherent can be obtained. This phase possesses a spin gap. We show that in this phase 
the order parameter for the BEC and the BCS superfluidity order parameter are identical, while charge density wave 
correlations present an exponential decay. We discuss the phase transitions induced by the detuning, the magnetic 
field and the repulsion. We also exhibit a solvable point where some correlation functions can be obtained exactly. 
In Sec. IIVI we consider the case where the number of atoms per site in the optical lattice is integer. We show that a 
phase transition to a Mott insulating state can be obtained in that case. However, there is no density wave order in 
this Mott state. Finally, in Sec. |Vj we discuss the applicability of our results to experiments. 



II. HAMILTONIANS AND THERMODYNAMICS 



A. Hamiltonians 

We consider a system of ID fermionic atoms with a Fano-Feshbach resonance. 0, I29I l80| This ID system can be 
obtained by trapping the fermions in a two dimensional or a three dimensional optical lattice. In the first case, the 
fermions are trapped into ID tubes, in the second case, a periodic potential is superimposed along the direction of 
the tubes. In the case in which the fermions are injected in a uniform potential, the Hamiltonian of the system reads: 
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Vbb(x - x')pb{x)pb{x') + Vff(x - x') 2J Pa(x)pa'(x') + 2V B f(x - x') 2J p a (x)pb(x') , (1) 

(7, a' cr 

where ipb annihilates a molecule, a fermion of spin cr, mp is the mass of the isolated fermionic atom, raj = 2mf 
the mass of the molecule, Vbb,Vbf,Vff are (respectively) the molecule-molecule, atom- molecule and atom-atom 
interactions. Since these interactions are short ranged, it is convenient to assume that they are of the form V a p(x) = 
g a p5(x). The term v is the detuning. Finally, the term A allows the transformation of a pair of fcrmions into a 
Fano-Feshbach molecule and the reverse process. This term can be viewed as a Josephson coupling |8lj between the 
order parameter of the BEC of the molecules, and the order parameter for the superfluidity of the fermions. As 
a result of the presence of this term, pairs of atoms are converted into molecules and vice-versa, as in a chemical 
reaction [Sfij. As a result of this, only the total number of atoms (paired and unpaired), M — 2Nb + Nf (where Nb is 
the number of molecules and Nf is the number of unpaired atoms) is a conserved quantity. 

In the case where atoms are injected in a periodic potential, V(x) — Vosiv? (iTx/d) it is convenient to introduce the 
Wannier orbitals[82j of this potential. In the single band approximation the Hamiltonian reads: [s^L 185ll8r3 | 

E( 6 V, + %+i) + u'j^M 2 +»Y, b h 

+AE( 6 i^T/id + fhfhbj) + V bf £ n b j{n m + n fdtl ), (2) 
j j 

where /j j(J annihilates a fermion of spin a on site j, nfj. a = fj a fj,a, bj creates a Fano-Feshbach molecule (boson) 
on the site j, and nbj — b^bj. The hopping integrals of the fermions and bosons are respectively t and t' . The 
quantity v is the detuning. The parameters U, U' and Vbf measure (respectively) the fermion- fermion, boson-boson, 
and fermion-boson repulsion. The case of hard core bosons corresponds to U' — > oo. The conversion of atoms into 
molecules is measured by the term A. Again, only the sum Af = 2Nb + Nf is conserved. We note that within the 
single band approximation, there should exist a hard core repulsion between the bosons. 




B. Thermodynamics of the boson-fermion model in the limit of A — > 

In this Section, we wish to study the behavior of the density of unpaired atoms pf and of the density of atoms 
paired in molecules pb as a function of the total density of atoms (pair and unpaired) ptot. in the limit of A — > 0+. In 
such a limit, the fermion-boson conversion does not affect the spectrum of the system compared to the case without 
fermion-boson conversion. However, it is imposing that only the total total number of atoms Af = 2Nb + Nf is 
conserved. Therefore, it this limit there is a single chemical potential p, and the partition function reads: 

Z x [p] = ^ mx ~ KNf+2Nb) \ (3) 

and: 

1 r)7 

N > + 2N * = W>^ (4) 

In the absence of fermion-boson conversion, A^f, and Nf would be separately conserved, and one would have a chemical 
potential pb for the molecules and pf for the atoms. The partition function of this hypothetical system would read: 

Z [pf,p b } = Tr{e-M H °-»f N f-^], (5) 

and thus: 

lira Z x [p] = Z [p,2p}. (6) 

A — >0+ 

If we further assume that Vbf = 0) we have H — Hf + Hb, where Hf is the Hamiltonian of the fermion subsystem and 
Hb is the Hamiltonian of the bosonic subsystem, and the partition function (jSJ) factorizes as Z [pf, pb] = ^ f[pf]Zb[pb], 
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where Z u [(x v ] = Tr^-^" - ^"] for v = f,b. Thus, in the limit A, Vbf -> 0, we obtain the following expression of 
the number of unpaired atoms Nf and the number of atoms paired in molecules Nf,. 



1 fdZ b \ 

We now use these equations Q and (JSJ to study the coexistence of bosons and fermions as the detuning v is varied. 
Two simple cases can be considered to illustrate this problem of coexistence. First, one can consider bosonic molecules 
with hard core repulsion and noninteracting fermionic atoms. In such a case, the thermodynamics of the gas of 
molecules is reduced to that of a system of spinless fermions by the Jordan- Wigner transformation)^ El) 113 > and the 
expression of the densities of unpaired atoms and molecules can be obtained in closed form. In this simple case, it is 
straightforward to show that for sufficiently negative detuning all atoms are paired into molecules, and for sufficiently 
positive detuning all the atoms remain unpaired. The case of intermediate detuning is more interesting as coexistence 
of unpaired atoms with atoms paired into molecules becomes possible. The physical origin of this coexistence is of 
course the molecule-molecule repulsion that makes the chemical potential of the gas of molecules increase with the 
density so that in a sufficiently dense gas of molecules, it becomes energetically favorable to create unpaired atoms. 
To show that the above result is not an artifact of having a hard core repulsion, we have also considered a slightly 
more realistic case of molecules with contact repulsion and non-interacting atoms. Although in that case we cannot 
anymore obtain closed form expressions of the density of molecules, we can still calculate numerically the density of 
molecules using the Lieb-Liniger solution 88] . We will see that having a finite repulsion between the molecules indeed 
does not eliminate the regime of of coexistence. 



1. The case of bosons with hard core repulsion 

In that case we assume that the boson-boson repulsion U' in the lattice case and gBB m the continuum case is 
going to infinity. Using the Jordan- Wigner transformation 87] , one shows that the partition function of these hard 
core bosons is equal to that of free spinless fermions thanks to the Jordan- Wigner transformation. For positive 
temperature, the number of unpaired atoms and the number of atoms paired into molecules read: 

Nt = 2L /— R , , (9) 

f dk 1 



For T — > 0, these equations reduce to: 



Pf = 
Pb = 



N f _ 2k F 

T ~ ~ : 

N^^kB_ 

L it ' 



fl = € F {kp) = , (11) 

where kp is the Fermi momentum of the atoms and kg is the Fermi momentum of the spinless fermions (i.e. the 
pseudo- Fermi momentum of the molecules). Up to now, we have not specified the dispersion of the atoms and of the 
molecules. In the lattice case, these dispersion are obtained from Eq. J5J) as ej (k) = —2t cos(k) and e&(fc) = —2t' cos(fc). 
A graphical solution of l|ll|l is shown on Fig. ^ for three different values of the chemical potential /x and v > 0. Three 
different regimes are obtained. In the first one, for = /x^, only unpaired atoms are present. In the second one for 
/.t = /is, unpaired atoms and molecules coexist. In the last one, for /i = /ic, all the available levels of unpaired atoms 
are filled, and the available levels for molecules are partially filled. As a result, the system behaves as if only molecules 
were present. This last phase is in fact a degenerate Tonks- Girardeau gas of molecules 0, In the intermediate 
regime, the fermions form a two-component Luttinger liauid|53.l55| and the bosons form a single component Luttinger 
liquid. 60] Similar calculations can be performed in the case of fermions and bosons in the continuum described by 
Eq. With free fermions and hard core bosons in the continuum Eqs. (|llfl become: 

kp 

2mp ' 
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1.2 
h D 



4m 



F 



v = 2/x, 



gPtot. = (fcF + fcs), (12) 

with p tot = 2/3^ + P_f the total density of atoms. Eliminating k F in Eq. 112|) . the problem is reduced to solving a 
second degree equation: 



3k B - 4npk B + n p - Am F v = 0. (13) 



The solutions of Eq. (|13|l arc: 



k F = \\/it 2 p 2 + l2m F v - ^p, 
3 6 

k B = 2 ^'^ V + 12m -^, (14) 
3 

and these solutions are physical when they yield both k F and k B positive. For v > 0, Eq. (|14fl yields k B > provided 
Ptot. > Ptot c = \\J m FV- When p < p[^ t c , the density of molecules is vanishing and ptot. = Pf- Above the critical 
density, atoms and molecules coexist, with densities given by Eq. I|12l) . At the critical density, the slope of k B versus 
p is discontinuous, being below the critical density and ? above the critical density. The Fermi wavevector k F also 
possesses a slope discontinuity at the critical density, the slope being zero above the critical density. The behavior of 
k F and k B as a function of the density is represented on Fig. 

For v < 0, Eq. (|T4fl yields k F > provided p > p^ t c — ^^Jm F \v\. When, p < p[^ t ^ c the density of unpaired atoms 
vanishes, and p — pb- Above the critical density, atoms and molecules coexist with densities given by Eq. (|12fl . As 
before, the slope of the curve k F versus p is discontinuous at the critical density, being zero below and 7r/3 above. 
The behavior of k F and k B as a function of the density for v < is represented on Fig- EJ 

The slope discontinuities in k B and k F have important consequences for the compressibility. Indeed, using Eq. I|12f) . 
it is easy to see that above the critical density, the chemical potential varies as 0(p— ptot.c) 2 - Since the compressibility 
X is defined as 1/x = p 2 |p , this implies that the compressibility of the system becomes infinite as the critical density 
is approached from above, signalling a first-order phase transition. Such first order transitions associated with the 
emptying of a band have been analyzed in the context of Luttinger liquid theory in Refs. j39, 90] . 

2. The case of bosons with finite repulsion 

We have seen in the previous Section that in the case of hard core repulsion between the molecules, both in the 
lattice case and in the continuum case, that having v < did not prevent the formation of unpaired atoms provided 
the total density of atoms was large enough. This was related with the increase of the chemical potential of bosons 
as a result of repulsion when the density was increased. In this section, we want to analyze a slightly more realistic 
case where the repulsion between bosons is finite and check that coexistence remains possible. In the lattice case, the 
problem is untractable by analytic methods and one needs to rely on numerical approaches . |9 ll l92l ] In the continuum 
case, however, it is well known that bosons with contact repulsions are exactly solvable by Bethe Ansatz techniques. jl3 
The density of molecules can therefore be obtained by solving a set of integral equations. |88l l93| They read: 

h 2 k 2 c f qo dq 

€{k) = ^ B - + "-» B+ n J_ go c 2 + (q-k) 2 eiq) > (15) 
90 P(q)dq 
-go c 2 + (k- qf 



2np(k) = l + 2cl -af^f-^ , (16) 



where: 



_ rn B g B B 



(17) 



/<?0 
p(q)dq, (18) 
-qa 

(]bb being the boson-boson interaction defined in Eq. The parameter q plays the role of a pseudo Fermi 

momentum. For q > qo, we have p(q) = 0. We also have e(±qo) = 0.|88| It is convenient to introduce dimensionless 
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FIG. 1: The different cases in the Bose-Fermi mixture with hardcore bosons when v > 0. For [i = fj,A, the Fermion band is 
partially filled, and the hardcore boson band is empty. For jj, = jj,c, the hardcore boson band is partially filled, and the fermion 
band is totally filled. In the case \x = ^b, both band are partially filled. In the rest of the paper we will only consider the latter 
case. 



variables 881: 



c c 
A = — ; 7 = — • 
qo PB 



h 2 2 

and rewrite k = qox, q — q^y, p(qox) = g{x), e(qox) — -^-^{x). The dimensionless integral equations read: 

2m{v — fis) A f 1 dy 



e(x) = x 2 + 



2ng(x) 



1 + 2A 



h 2 q 2 7T J_ x \ 2 + {x- y) 2 

g{y)dy 



_ x \i + (x-y) 2 ' 

Using e(±go) = one has the following integral equation for e(x): 



e(x) = x 2 - 1 + 



A 



dye{y) 



X 2 + (x ~ y) 2 A 2 + (1 — y)' 2 



Once this equation has been solved, the chemical potential of the bosons is obtained by: 



Ms = v 



2m B 



1 



7T J_x X 2 + (X - 1) : 



;e(x)dx 



(19) 

(20) 
(21) 

(22) 

(23) 
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7ip 

FIG. 2: The behavior of &f and fes for positive detuning v > as a function of the total density p. For low densities, only 
atoms are present (fcs = 0). At higher densities such that 7rp > 2^mpu, a nonzero density of molecules appear. At the critical 
density, the slopes of Uf and fcs versus p are discontinuous. On the figure we have taken vtlfv = 1. 




4 5 6 7 8 9 

FIG. 3: The behavior of ItF and fcs for negative detuning < 0. For low densities, only the molecules are present (fej? = 0). 
For 7rp > m^M, molecules coexist with atoms. At the critical density, the slopes of Uf and ks versus p are discontinuous. 
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Knowing ps gives immediately (If = Ms/2. From pp one finds Uf — \f2rappp and pp = Ikp/ir. Finally, using the 
definition of the total density p = 2ps + Pf one can map the molecule density and the free atom density as a function 
of the total density of atoms. The resulting equation of state can be written in terms of dimensionless parameters as: 

h 2 p B _ <p ( ^ v 

niBgiD \ m B giD ' rn F gl D 

The behavior of the boson density ps and fermion density pf as a function of total density p can be understood 
in qualitative terms. Let us first discuss the case of negative detuning. For sufficiently low densities, only bosons are 
present. However, in that regime, the boson-boson repulsion is strong, and the boson chemical potential is increasing 
with the boson density. As a result, when the density exceeds a critical density p Cl the fermion chemical potential 
becomes positive, and the density of fermions becomes non-zero. The appearance of fermions is causing a cusp in 
the boson density plotted versus the total density. When the density of particles becomes higher, the boson-boson 
interaction becomes weaker, and the boson chemical potential barely increases with the density. As a result, the 
fermion density becomes almost independent of the total density. In the case of positive detuning, for low density, 
only fermions are present. Again, the increase of fermion density results in an increase of chemical potential and 
above a certain threshold in fermion density, bosons start to appear, creating a cusp in the dependence of the fermion 
density upon the total density. At large density, the detuning becomes irrelevant, and the fermion density barely 
increases with the total density. 

To illustrate this behavior, we have solved numerically the integral equations 120(1. and calculated the resulting 
fermion and boson densities. A plot of the density of bosons as well as the density of fermions is shown on Fig. 0] for 
v > and on Fig. [S] for v < 0. The slope discontinuities at the critical density remain visible. Obviously, this implies 
that the divergence of the compressibility is still present when the repulsion between the molecule is not infinite. 
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FIG. 4: The density of molecules ps and unpaired atoms pf as a function of the total density p in the case of a repulsion 
c = 100 between the bosons and for positive detuning v — 0.1. At large density, the fermion density is increasing more slowly 
than the boson density. Inset: the behavior of the boson and fermion densities near the origin. Note the cusp in the fermion 
density as the boson density becomes nonzero as in the c = oo case. 

We have thus seen that generally we should expect a coexistence of fermionic atoms and bosonic molecules as soon 
as repulsion between the molecules is sufficiently strong. Moreover, the repulsion between the molecules results in a 
finite velocity for sound excitations in the molecule Bose gas. As a result, we can expect that the gas of molecules 
will behave as a Luttinger liquid. Till now however, we have assumed that the term converting atoms into molecules 
was sufficiently small not to affect significantly the spectrum of the system. In the following, we will treat the effect 
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FIG. 5: The density of molecules ps and unpaired atoms pf as a function of the total density p in the case of a repulsion 
c = 100 between the bosons and for negative detuning v — —0.1. At large density, the fermion density is increasing more slowly 
than the boson density. Inset: the behavior of the boson and fermion densities near the origin. Note the cusp in the boson 
density as the fermion density becomes nonzero as in the c = oo case. 

of a small but not infinitesimal conversion term in Eqs. J5J and Q using bosonization techniques. We will show that 
this term can lead to phase coherence between the atoms and the molecules, and we will discuss the properties of the 
phase in which such coherence is observed. 



In this Section, we consider the case discussed in Sec. ITT1 where neither the density of molecules nor the density of 
atoms vanishes. As discussed in Sec. [HJ this requires a sufficiently large initial density of atoms. As there is both a 
non-zero density of atoms and of molecules, they both form Luttinger liquids |54l l5ol l6Cl | . These Luttinger liquids are 
coupled by the repulsion between atoms and molecules Vbf and via the conversion term or Josephson coupling A. 
To describe these coupled Luttinger liquids, we apply bosonization [5^ to the Hamiltonians QJ. For the sake of 
dcfiniteness, we discuss the bosonization procedure in details only in the case of the continuum Hamiltonian Jy|. For 
the lattice Hamiltonian J5J, the steps to follow are identical provided the system is not at a commensurate filling. At 
commensurate filling, umklapp terms must be added to the bosonized Hamiltonian and can result in Mott phases |53|. 
This case is treated in Sec. IIVI 

To derive the bosonized Hamiltonian describing the low-energy spectrum of the Hamiltonian we need first to 
consider the bosonized description of the system when all atom-molecule interactions are turned off. For A = 0, Vbf = 
0, both Nf and Nb are conserved and the bosonized Hamiltonian equivalent to or 0J is given by: 



III. PHASE DIAGRAM AND CORRELATION FUNCTIONS 



A. 



Derivation of the bosonized Hamiltonian 



H 



Hb + Hp + H a 




u b K b (7rUb) 2 + -±(d x <k) 2 
u p K p {ttTI p ) 2 + ^f{d x 4> P ) 2 




p 
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Ha = UaK^U^ 2 + ^-(d^) 2 -JfHa JdxaxVfy, (25) 

where [^(x). IL/ (a/)l = i<5(x — x')8 v y , (i/, i/ = b, a, p). In the context of cold atoms, the Hamiltonian (|25|l have been 
discussed in |54L I55L mOl] . The parameters K a , the Luttinger exponent, and u p ,Ua-, the charge and spin velocities, are 
known functions of the interactions [toL l94l l95l | . with K p = 1 in the non-interacting case, gi± is a marginally irrelevant 
interaction, and at the fixed point of the RG flow K* — 1. For the bosonic system, the parameters Uh, Kb can be 
obtained from numerical calculations 92] in the lattice case or from the solution of the Lieb-Liniger model 88] in the 
continuum case. In the case of non- interacting bosons Kb — > oo and in the case of hard core bosons Kb = !.|59ll6lL l62]] 
An important property of the parameters Kb and K p is that they decrease as (respectively) the boson-boson and 
fermion-fermion interaction become more repulsive. The bosonized Hamiltonian (|25[1 is also valid in the lattice case 
(0) provided that both Nf and Nb do not correspond to any commensurate filling. 

The fermion operators can be expressed as functions of the bosonic fields appearing in i(23)l as |53l ]: 

Ms) = Yl e irkFna Mx = na) (26) 

r=± 

(x) = = , (27) 

where the index r = ± indicates the right /left movers, a is a cutoff equal to the lattice spacing in the case of the 
model Eq. (0). Similarly, the boson operators are expressed as[53|: 

-JL = ^ b (x = na) (28) 
V a 

^b(x) = -^L= [l + A cos(20 6 - 2k B x)\ . (29) 
V27ra 

In Eqs. (|26|l - 128|l . we have introduced the dual fields[5i| v {x) — ir f x TL v (x')dx' (v = p,a,b), Uf = nNf/2L, and 
ks = nNb/L where L is the length of the system. The fermion density is given bv|53|: 



E "<■"■■'■ , , . cos(2fc F x - V20 p ) r- 
= pf(x = na) — o x q>p-\ — cosv20 (T , (30) 
a tt na 

a 

and the boson density bv|53|]: 



a 7r 
The detuning term in is thus expressed as: 



n b , n f , 1„ , , cos(2k B x - 20b) 

Pb(x) = --d x 4>b H • (31) 



7T 



H detuning = / dxd x (j>b (32) 



We now turn on a small A and a small Vbf- The effect of a small Vbf on a boson- fermion mixture has been 
investigated previously [§E H3- The forward scattering contribution is: 



V BF V2 



d x 4>bd x 4> p , (33) 



and as discussed in [96j, it can give rise to a phase separation between bosons and fermions if it is too repulsive. 
Otherwise, it only leads to a renormalization of the Luttinger exponents. The atom molecule repulsion term also gives 
a backscattering contribution: 

2V BF 



{2tto)' 



J dxcos(2<j>b — V2~4> p — 2(kp — k B )x) cosv^cr, (34) 

however in the general case, kp ^ ks this contribution is vanishing. In the special case of k B = kp, the backscattering 
can result in the formation of a charge density wave. This effect will be discussed in Sec. IIII B 31 The contribution of 
the A term is more interesting. |76l 1771] Using Eqs. $ty - $2fy . we find that the most relevant contribution reads: 



TT 2A 
tl b f 



V27r 3 a 3 



(6 b -V28 p ) cos V2(t) a (35) 



In the next section, we will see that this term gives rise to a phase with atom-molecule coherence when the repulsion 
is not too strong. 
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B. Phase diagram 



phase with atom-molecule coherence 



The effect of the term (|35f) on the phase diagram can be studied by renormalization group techniques!^ • A detailed 
study of the renormalization group equations has been published in [76} . Here, we present a simplified analysis, which 
is sufficient to predict the phases that can be obtained in our system. The scaling dimension of the boson-fermion 
coupling term l|35[l is: j^- + + \K n . For small A it is reasonable to replace K a with its fixed point value K* = 1. 



Therefore, the RG equation for the dimensionless coupling A 



At. 1 ' 2 



(where u is one of the velocities u p , u ai u b ) reads: 



dX 



1 



2K„ 



(36) 



where t is related to the renormalized cutoff a(£) — ae i . We thus see that for — h < 3/2, this interaction is 
relevant. Since for hardcore bosons [6lt l62j Kb = 1 and for non-interacting bosons Kb = oo, while for free fermions 
K p = 1 and in the lattice case for U = oo one has K p = 1/2(9^, we see that the inequality is satisfied unless there 
are very strongly repulsive interactions both in the boson system and in the fermion s ystem . When this inequality 
is not satisfied, for instance in the case of fermions with nearest-neighbo r rep ulsion[99lli"on|. in which one can have 
1/4 < K p < 1/2 and hardcore bosons with nearest neighbor repulsion [H?! Il0lj | . in which one can have K b = 1 /2, the 
atoms and the molecules decouple. This case is analogous to that of the mixture of bosons and fermions |96l l97j]and 
charge density waves can be formed if ks and Uf are commensurate. The phase transition betwee n thi s decoupled 
phase and the coupled phase belongs to the Berezinskii-Kosterlitz-Thouless (BKT) universality class. [l02j As pointed 
out in |76|. in the decoupled phase, the effective interaction between the fermions can be attractive. In that case, a 
spin gap is formed [5^, and the fermions are in a Luther- Emery liquid state with gapless density excitations. Let 
us consider the coupled phase in more details. The relevance of the interaction (|35|l leads to the locking of i.e. it 
results in the formation of a spin gap. To understand the effect of the term cos(#b — \/29 p ), it is better to perform a 
rotation: 



(37) 



and the same transformation for the (j) v . This transformation preserves the canonical commutation relations between 
cj)± and II±. Under this transformation, Hi, + H p becomes: 




Hh + H n 



dx 
2^ 



E 



K u 



f dx 

+ J —[g 1 {^ + ){^)+g 2 d x <j 3+ d x 



(38) 



where: 



u+K+ 
U-K- 

.91 

u+ 
K+ 

U- 

k~Z 

.92 



u b K b 
u b K b 



^UpKp, 
2 

^UpKp, 



^-{u b K b - UpKp), 



2u b 



3K b 3K p 3vr 
u b 2u p 
3K b 3K P 



3^' 



V", 



while, H b f defined in (|35|l becomes: 



Hi 



VS U b , 
u„K n , 

3 ^Kl, P P 7T 



• / dxcos \/39_ cosV2</v. 



(39) 
(40) 
(41) 

(42) 
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After the rotation, we see that when A is relevant, the field 6L is also locked, but (f> + remains gapless. Since the field 
6>_ is the difference of the superfluid phase of the atoms and the one of the molecules, this means that when A becomes 
relevant, unpaired atoms and molecules share the same superfluid phase i.e. they become coherent. The gap induced 
by the term A can be estimated from the renormalization group equation l|3()|) ■ Under the renormalization group, the 
dimensionless parameter X(£) grows until it becomes of order one at a scale t — I* where the perturbative approach 
breaks down. Beyond the scale £*, the fields 6L and <j) a behave as classical fields. Therefore, the associated energy 
scale u/{ae l ) is the scale of the gap. From this argument, we obtain that the gap behaves as: 



A — — 



Xa 1 ' 2 



— r 



— r~ 

4A-,, 



(43) 



The gapful excitations have a dispersion e(k) = J (uk) 2 + A 2 and are kinks and antikinks of the fields 6L and ^. [lOS 
More precisely, since a kink must interpolate between degenerate classical ground states of the potential l|42[) . we find 
that when a kink is present 6>_(+oo) — 0-{— oo) = ±7r/\/3 and </><j(+oo) — cf> a (~oo) — n/\/2. This indicates that a 
kink is carrying a spin 1/2, and is making the phase 9b of the bosons jump by 7r/3 and the phase of the superfluid 
order parameter V29 p of the fermions jump by — 27r/3. Since the current of bosons is jb = UbKbirllb = UbKbd x 9b and 
the current of fermions is jp — V2u p K p ttII p = \/2u p K p d x 9 p , this indicates that counterpropagating supercurrents of 
atoms and molecules exist in the vicinity of the kinks. Therefore, we can view the kinks and antikinks are composite 
objects formed of vortices bound with a spin 1/2. We note that the kinks and antikinks may not exhaust all the 
possible gapful excitations of the system. In particular, bound states of kinks and antikinks, known as breathers may 
also be present [l03j . However, these gapful excitations present a larger gap than the single kinks. Let us now turn to 
the gapless field <p+. This field has a simple physical interpretation. Considering the integral 



dxd x 4> + 



(x/24 + ^) = ^=, 



V6' 



(44) 



showing that <f> + (oo) — 4>+(— oo) measures the total number of particles in the system N . Thus (II + , < 
total density excitations of the system. 

The resulting low-energy Hamiltonian describing the gapless total density modes reads: 



.) describe the 



"-='§ 



K + 



(45) 



where u*^_,K^_ denote renormalized values of u +1 K + . This renormalization is caused by the residual interactions 
between gapless modes and the gapped modes measured by gi , g 2 in Eq. I|38|) . Since + measures the total density, 
the Hamiltonian 145|) describes the propagation of sound modes in the ID fluid with dispersion uj(k) = u\k\. We note 
than in Ref. [7^. Il04| . dispersion relation similar to ours were derived for the sound modes and the superfluid phase 
difference modes using different methods. 



2. effect of the detuning and applied magnetic field 

Having understood the nature of the ground state and the low excited states when A is relevant we turn to the 
effect of the detuning term. Eqs. (|32|l and (|37|l show that the detuning term can be expressed as a function of 4>+, </>- 
as: 




This shows that the detuning does not affect the boson-fermion coupling (|35|l since it can be eliminated from 

the Hamiltonian by a canonical transformation tfi± — > <ft± + X±x, where A + = arid A_ = v ~^- For a fixed 

total density, changing the detuning only modifies the wavevectors ks and kp. As discussed extensively in Sec. Ill Bl 
for a detuning sufficiently large in absolute value, only molecules or only atoms are present, and near the critical 
value of the detuning, the compressibility of the system is divergent. We therefore conclude that in one-dimension, 
the crossover from the Bose condensation to the superfluid state, as the detuning is varied, is the results of the 
band-filling transitions at which either the density of the atoms or the molecu les g oes to zero. At such band filling 
transitions, v p a — > (respectively Vb — > 0) and bosonization breaks down |89ll9Ctll05j . The two band filling transitions 
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are represented on Fig.0 The cases at the extreme left and the extreme right of the phase diagram have been analyzed 
in [XJ], where it was shown that in the case of a broad Fano-Feshbach resonance, the zone of coexistence was very 
narrow. In the narrow Fano-Feshbach resonance case we are investigating, the zone of coexistence can be quite 
important. 

Band-filling transitions 



No fermions 'Coexistence No bosons 

1 1 ^ 

v 

FIG. 6: The band filling transitions as a function of the detuning 

Application of a magnetic field can also induce some phase transitions. The interaction with the magnetic field 
reads: 

h 

7lV2~, 

The effect of the magnetic field is to lower the gap for the creation of kink excitations (remember that they carry 
a spin 1/2). As a re sult, when it becom es larger than the gap, the magnetic field induces a commensurate incom- 
mensurate transition |l06l Il07l Il08l Il09f that destroys the coherence between atoms and molecules and gives back 
decoupled Luttinger liauids |ll(3| . In that regime, the behavior of the system is described by the models of Ref. |9^, In- 
commensurate-incommensurate transitions have been already discussed in the context of cold atoms in |lll| . In the 
problem we are consideri ng, however , since two field s are becoming gapless at the same time, and <j) a , there are 
some differences [Til EH with the standard case [m|, in particular the exponents at the transition are non-universal. 

To conclude this section, we notice that we have found three types of phase transitions in the system we are 
considering. We can have Kosterlitz-Thouless phase transitions as a function of interactions, where we go from a 
phase with locked superfluid phases between the bosons and the fermions at weak repulsion to a phase with decoupled 
bosons and fermions at strong repulsion. We can have band-filling transitions as a function of the detuning between 
the phase in which atoms and molecule coexist and phases where only atoms or only molecules are present. Finally, 
we can have commensurate-incommensurate transitions as a function of the strength of the magnetic field. In the 
following section, we discuss the correlation functions of superfluid and charge density wave order parameters in the 
phase in which molecules and atoms coexist with their relative superfluid phase 9- locked. 



H magn = 7^ / dxd x (j)<j, (47) 



3. Quantum Ising phase transition for kf 



In the case of kp = ks, the backscattering term (|34|l is non- vanishing. This term induces a mutual locking of the 
densities of the bosons and the fermions 96] and favors charge density wave fluctuations. This term is competing 
with the Josephson term l|35|l which tends to reduce density wave fluctuations. For fei? = ks the relevant part of the 
Hamiltonian given by the combination of the terms (|34|l and l|35f) reads: 



H 



Joscphon+CDW L. 



ock. / dx 



J= cos{9 b - s/26 p ) + cos(20 fc - V2> p 



Using a transformation </>{, = 0J./V2, 9b — 9bV2, and introducing the linear combinations 

, <f>b + <f> P 

± 4>b- <t>p 



»2 = 



V2 



and similar combinations for the dual fields, we can rewrite the interaction term l|48|l as: 



H. 



Joscphon+CDW Lock. 



— dx 



2A 2V BF 

cos 29 2 + — — — cos 202 



y/2ir 3 a 3 



(27ra) 



cos 



V20 o 



(48) 

(49) 
(50) 

(51) 
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From this Hamiltonian, it is immediate to see that a quantum Ising phase transitio n occurs between the density wave 
phase 4>2 and the superfluid phase 8 2 at a critical point A c = -^y^. jrLlL Ill5l Ill6| Indeed, the field 4> a being locked, 

we can repla ce cos \/2(t> a by its expectation value in Eq. I]5ip. and rewrite (|51|l as a free massive Majorana fermions 
Hamiltonian. |ll4 . Ill5l lrTa| At the point X n . the mass of one of these Majorana fermions vanishes giving a quantum 
critical point in the Ising universality class |ll7| . On one side of the transition, when A > A c , the system is in the 
superfluid state discussed in Sec. IIII Bl on the other side A < A c , the charge density wave state discussed in [23 is 
recovered. 



Correlation functions 



In order to better characterize the phase in which A is relevant, we need to study the correlation function of the 
superfluid and the charge density wave operators. Let us begin by characterizing the superfluid order parameters. 
First, let us consider the order parameter for BEC of the molecules. As a result of the locking of the fields 6L and 
(j> a , the boson operator Eq. (|28|l becomes at low energy: 



*s(x) 



3°+ 



(52) 



An order of magnitude of (e ' can be obtained from a scaling argument similar to the one giving the gap. Since 

. e_ 

the scaling dimension of the field e * "75 is 1/12K-, and since the only lengthscale in the problem is e l , we must have 
.»_ 

(e~ 1 ^) ~ e~ e / 12K - ^ [Xa 1 / 2 /u) 1 / 12K - . Similarly, the order parameter for s-wave superconductivity of the atoms 
Oss = Yla^r^-r-a becomes: 



o 



ss 



£ V3 



(53) 



thus indicating that the order parameters of the BEC and the BCS superfluidity have become identical in the low 
energy limit |3lL l32j|. This is the signature of the coherence between the atom and the molecular superfluids. The 
boson correlator behaves as: 



OMar,T)tf B (0,0) 



a 



x 2 + (ut)' 



(54) 



As a result, the molecule momentum distribution becomes ns(fc) ~ |fc| 1 /( 3 - K +) -:1 . One can see that the tendency 
towards superfluidity is strongly enhanced since the divergence of ns{k) for k — > is increased by the coherence be- 
tween the molecules a nd the ato ms. This boson momentum distribution can, in principle be measured in a condensate 
expansion experiment |ll8l Ill9| . 

Having seen that superfluidity is enhanced in the system, with BEC and BCS order parameters becoming identical, 
let us turn to the density wave order parameters. These order parameters are simply the staggered components of 
the atom and molecule density in Eqs. (|30l) — (|3 1|) . In terms of <p±, the staggered component of the molecule density is 
reexpressed as: 



P2k B ,b( X ) 



cos 




2 ^= + ^-=(j) + - 2k b x 



and the staggered component of the fermion density as 

p2k F ,f(x) ~ cos 




- 2k FX 



(55) 



(56) 



where we have taken into account the long range ordering of </> CT . We see that the correlations of both p2k B ,b( x ) an d 
P2kpj( x ) decay exponentially due to the presence of the disorder field 0_ dual to 6>_. In more physical terms, the 
exponential decay of the density-wave correlations in the system results from the constant conversion of molecules 
into atoms and the reciprocal process which prevents the buildup of a well defined atomic Fermi-surface or molecule 
pseudo-Fermi surface. The exponential decay of the density wave correlations in this system must be be contrasted 
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with the power-law decay of these correlations in a system with only bosons or in a system of fermions with attractive 
interactions. |53j In fact, if we consider that our new particles are created by the operator ipf, ~ V'tV 1 !; we can derive 
an expression of the density operators of these new particles b y con sidering the product ^bV'T^l- Using the Haldane 
expansion of the boson creation and annihilation operators |59l Tl20j . we can write this product as: 



cos(2m</> — 2mkBX 

.m=0 
-V2B P ) 



x e 



i\/28 p 



cos V2(f> a + cos(V2(f)p — 2kFx) 



cos 



V2<t>a) cos(V6<ft+ — 2{k F + k B )x), 



(57) 



where {k F + ks) — Tt{2N[, + Nf)/2L = nppairs can be interpreted as the pseudo Fermi wavevector of composite 
bosons. A scaling argument shows that the prefactor in the expression varies as a power of A. As a result, when there 
is coherence between atoms and molecule, power-law correlations appear in the density-density correlator near the 
wavevector 2k F + 2ks and the intensity of these correlations is proportional to the | (er l ( eb ~^ e ^ cos \/2(j)a) \ 2 ■ The 
resulting behavior of the Fourier transform of the density-density correlator is represented on Fig. 



Q=2k 



B 



Q=2k F 



Q=2k B + 2k F 




O 

FIG. 7: Fourier transform of the static density density correlations. In the decoupled phase (dashed line), two peaks are 
obtained at twice the Fermi wavevector of the unpaired atoms and at twice the pseudo-Fermi wavevector of the molecules. In 
the coupled phase (solid line), the peaks are replaced by maxima at Q = 2k F and Q = Iks- A new peak at Q = 2(k F + fcs) is 
obtained as a result of Boson-Fermion coherence. 



Another interesting consequence of the existence of atom/molecule coherence is the possibility of having non- 
vanishing cross-correla tion s of the atom and the molecule density. In the three-dimensional case such cross correlations 
have been studied in |l2l| . If we first consider cross-correlations (T T p2k B ,b(x,T)p2k F ,f(0,0)} we notice that due to 
the presence of different exponentials of in Eqs. (|55|l - 1)56(1. this correlator vanishes exactly. Therefore, no cross 
correlation exists between the staggered densities. However, if we consider the cross correlations of the uniform 
densities, we note that since they can all be expressed as functions of d x (j) + ,d x (j)-, such cross correlations will be 
non- vanishing. More precisely, since: 



PF 



V2_ 
7T V3 



d x (j)- 



(58) 
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PB 



V2_ 



ir V3 



(59) 



at low energy we have: pf ~ Pb ~ ~~0?P X 



D. The Luther Emery point 

In this Section, we will obtain detailed expressions for these correlation functions at a special exactly solvable point of 
the parameter space. At this point, the kinks of the fields 4>a and 0- become free ma ssive fermions. This property , and 
the equivalence of free massive fermions in ID with the 2D non- critical Ising model [IH EH EM EH uM, EH Ell 
allows one to find exactly the correlation functions. 



1. mapping on free fermions 

As we have seen, after the rotation 1)3 7|l. if we neglect the interaction terms of the form n+II_ or d x 4>+d x 
Hamiltonian of the massive modes </>_, <j) a can be rewritten as: 



the 



H 



dx 
2^ 



+ 



dx 
2^ 
A 



«_iC_(7rn_; 



-{d x <t>-f 



•\/27r 3 a 3 

where K* = 1. When the Luttinger exponent is K— 



dx COS VS9- COS \f2(j> a , 

3/2, it is convenient to introduce the fields: 



(60) 

(61) 
(62) 



and rewrite the Hamiltonian H6()|l as: 



H 



J J K(»n.)' 



dx 
2^ 
A 



[ u _(^n) 2 + «_(a^) 2 ] 

/dx cos \/20 cos \/2<j) a . 



If we neglect the velocity difference, i.e. assume that u* = u_ = it, and introduce the pseudofermion fields: 

e i[(e-r^)+cr(0„-r0„)] 
*r .„ = — 



\/2ira 

we see immediately that the Hamiltonian l|63[) is the bosonized form of the following free fcrmion Hamiltonian: 

A 



H 



E 



dx 



-iu r*J: i<r fl B *r,<7 + -j== 



(63) 



(64) 



(65) 



As a result, for the special value of K- = 3/2, the excitations can be described as massive free fermions with 
dispersion e(k) = \J (uk) 2 + m 2 , where m — -y . This is known as Luther-Emery solution. |58l Il29^ One can see 
that the fermions carry a spin 1/2 and a jump of the phase equal to Therefore they can be identified to the 

kinks obtained in the semiclassical treatment of Sec. IIII 51 Also, making all velocities equal and Vbf = in l|39|l . 

i 

we find the relation 3/K- = 1/K b + 2/K p and thus the gap given by the RG varies as A ~ ^ ■ For 

K- =3/2 this expression reduces to the one given by the fermion mapping. 
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2. Correlation functions 



To calculate the correlation functions it is convenient to introduce the fields: 



(66) 
(67) 



And reexpress operators by: 



= e i(e T +ei) 



-.T^+a 



This leads us to the following expression of the fermion density p2k F , /(%)'■ 



(68) 



(69) 



which enables us to find exactly its correlations. We introduce u = u* to simplify the notations. The Green's functions 
of the fermions read: 



G++(x 1 t) = - 
G-(x-t) = - 



m 


T + 


U 


2nu 




+ 


m 


T — 


U 


2ttu 




+ 



^=Ki m\ t 2 + 



G- + (x, t) = G+ _ (aj>r) = _^o(m^ + (2) 3 ) 



(70) 
(71) 
(72) 



Using (|69(l and Wick's theorem we find that in real space the density density correlations read 

(T T p2k F j(x,T)p-2k F ,f(0-0)) = 2 



V 2ttu ) \x 2 + (ut) 2 



Ki \m\ t 



K(\m\r 2 



(73) 



These correlation functions decay exponentially, with a correlation length u/m = £. Note that expression 17311 is 
exact. 

On the other side the boson density is given by: 



P2k B ,b 



(74) 



To calculate the correlation fu nctions in this case , we can use the equivalence of the Dirac fermions in (1+1)D with 
the 2D non critical Ising modelp l IrM ITM H2H fl27l | to express the boson fields in terms of the order and disorder 
parameters of two non-critical Ising models, 0x2,^1,2 respectively, by: 

(75) 
(76) 
(77) 
(78) 

We find that: 



(T T P2k B ,b(x, T)p 2 k B ,6(0, 0)) 



COS0 CT 


= <Ji p 2 


sin <3 a 


= Ml 02 


cos $ CT 


= °~lO~2 


sin$ CT 


= P-lP-2 


a 2 


) 3 ^ 


+ {ut) 2 



4(a(x,T)a(0,0)) 2 (Mx,T)/i(0,0)r 



(79) 



18 



where we have used (ai,2 {x, t)(7i )2 (0, 0)) = (cr(x, r)a(Q, 0)) and similarly for /i. In the boso nic c ase, the mapping 
on the 2D Ising model allows to calculate the correlation functions using the results of Ref. Il3f|. w here an exact 
expression of the correlation functions of the Ising model in terms of Painleve III transcendants [l31| was derived. 
In fact, since we are interested in the low-energy, long-distance properties of the system, it is enough to replace 
the Painleve transcendants with an approximate expression in terms of modified Bessel functions. The resulting 
approximate expression is: 



2K + 

(T T p 2fcB , 6 (s,r) P2feB , b (0,0)) ~ ( ^ 2+ ^ r)2 ) 3 Kl L^ T 2+^) 2 j. (80) 

Knowing the correlation function in Matsubara space allows us to obtain them in the reciprocal space via Fourier 
transforms. The Fourier transform of the density-density response functions l|73[) and (|80|l can be obtained in Mat- 
subara space from integrals derived in the Appendix lAl 

We find that the bosonic structure factor is: 



x£(±2* fl + g>w ) = -(_) (-) ^ 3 __ 2J J 3^(1-^1-^1-^^-^!;- 



3 



2 



4m 2 



where r(x) is the Gamma function and 3-F2C z ) is a generalized hypergeometric function. |132j Of course, 

since Eq. I|8(J|I is a long distance approximation, the expression l|81|) is also approximate. More precisely, the exact 
expression possesses thresholds at higher frequencies associated with the excitation of more than one pair of kinks in 
the intermediate state. However, the expression H81fl is exact as long as w is below the lowest of these thresholds. For 
the fermions, the expression of the structure factor is exact and reads: 



w 1 , ~i \ 1 rma\ ■ 



K + 
3 










r 


(i 





tf+ JT+ 3 X + c^ 2 + ( M<z ) : 
3^2 [ l — , l — j l — i 7; —■> x ' _ " 



6 ' 6 ' 6 ' 2 6 ' ' 4m 2 



r(2- ^)r(i-g*)r(-^) / k+ _ k+ _k+, 3 _ jc+ c 2 + ( Mg ) 2 

r (3_^\ 32 V 6' 6' 6 '2 6' ' 4m 2 lq 



The response functions are then obtained by the subs titution ioj — > + iO. Since the generalized hypergeometric 
functions p+ ii ? p (. z) are analytic for \z\ < 1 |133| |. the imaginary part of the response fun ctions is v anishing for 
lu < 2m. For u) > 2m, the behavior of the imaginar"rt is obtamed from a theorem quoted in EH US- According 
to the theorem, 



T m\'''m P ) p+lFp{ai ■ ■ ■ ap+i;bl ■ ■ ■ bp;z) = £ 5m(0)(1 " zr + (1 " z) " P £ ffm(s?,)(1 " z)m ' (83) 



provided that: 



v p+i 



i=l i=l 

is not an integer. Therefore, if s p < 0, the generalized hypergeometric function possesses power-law divergence as 
z — ► 1. If < s p < 1, it has a cusp singularity. 
In our case, for the fermion problem, we have 

r F_, , 3 K + f, K + \ _ K + 1 

and in the boson problem 



<- 1 + 5-f-H 1 -fJ-f-5 (85) 



19 



Therefore, for K + small, both the fermion and the boson density-density response functions show power law singulari- 
ties for uj — * 2m. For K + = 3/8, the singularity in the boson density-density response is replaced by a cusp. This cusp 
disappears when K + = 9/8. For K + = 3/2 the singularity in the fermion density density correlator also disappears. 
The behavior of the imaginary part of the boson density-density correlation function is shown in Fig|S| 

An ex act expression of the imaginary part of the 3-F2 function can be deduced from the calculations in |l34j . Indeed, 
in |134| . it was found that: 

t? I 7 j \ tp t 1 7 \ . r(6i)r(oi + a 2 + q 3 - bi - 6 2 ) tl , , , , . 

3 F 2 (ai, a 2 , a 3 ; 61, 62; z) = F R (ai,a 2 , a 3 ; 61, 6 2 ; z) H — — — — — - — - t(ai, a 2 , a 3 ; b 1 ,b 2 ; z), (87) 

r(ai)r(a 2 )r(a 3 ) 

where Fr is defined by a series that converges absolutely for Re(z) > 1/2 and thus has no cut along [1, +00] , and f; is 
singular along [1, +00]. £ can be expressed in terms of a higher hypergeometric function of two variables, the Appell 
function F 3 , as: 

£(01,03,03; 61,62;*) = Z ^-b,+l^_ z f,+b,-a t -a 2 -a 3 

xF 3 (b 1 — oi, 1 - 02, 62 — oi, 1 — a 3 , 61 + 62 — oi — 02 — 03 + 1, 1 — 1/z, 1— z). (88) 

Since only £ has a cut along [l,oo], for u> 2 > (vq) 2 + 4m 2 , z > 1, the imaginary part can be expressed as a function 
of F3 only. This is particularly useful in the case of the fermion density correlators, because the expression H82JI is 
exact. In the case of the bosons, thresholds associated with the excitation of a larger number of Majorana fermions 
will appear a t energies ~ 4m . The imaginary parts of correlation functions Eq. H81II - (|82|) can be measured by Bragg 
spectroscopy [Hi Ellini- In FigElwe plot the imaginary part of density correlation functions for the fermionic 
system with K + = 1/4 and K + = 1/2 as a function of frequency. In the first case, we obtain a divergence of the 
density-density response near the threshold, whereas in the second case, only a cusp is obtained. 



3. spectral functions of the fermions 



We also wish to calculate the spectral functions of the original fermions Vv,<r (not the pseudofermions ^r,a)- To 
obtain these spectral functions, we express the operators ip r ,a as a function of the fields <p+, and their dual fields. 
We find: 



1 _!_ 



'2ira 

1 =e 7f( e ++^+) e 4 e -+l*<'] e -i*- 



'2na 



Therefore, the Green's functions of the original fermions factorize as: 

- (T t ^ + .„(x,t)^„(0,0)) = G + (x,r)GZ 1 f /6 (x,T)G° a 1 / 6 (x,r), 

where: 



G+(x,t) = -(T T eV5 



(6+-<t>+)(x,T) e -^(0 + -0+)(O : O)> 



UT — IX \X 2 + (ut)' 2 



(89) 
(90) 

(91) 
(92) 



GI^ 5/6 (x,r) = _( Tr e4- e -+l*-](^-) e 4^-|*^(o.o)), (93) 

and: 

G^lTix.T) = -(T T e-^ x ^e^^). (94) 



The correlator in Eq. I|94l) satisfies differential equations that were derived in |l39j| . However, since the fields & a are 
long range ordered, we can simply replace the terms e by their expectation value (e This approximation 

only affects the behavior of the fermion correlator at high energy. 

Therefore, we are left with l|9"5|) to evaluate. To do this, we can use exact results derived in |l40Ul4llll42fl to obtain 
a form-factor expansion of Eq. (|93(l . The first term of the Form factor expansion yields: 
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FIG. 8: The imaginary part of the density-density correlation function for the fermionic system with K+ — 1/4, 1/2. 
Writing: 



ut = pcoscp (96) 
x = psimp (97) 



We finally obtain that: 



To [K 4 



<Z(*,r)~e*^j + K^p). (98) 

The Fourier transform of G(x, t) is given by a Weber-Schaefheitlin integral pTl [Til. [Til Hil Il45| . In final form, 
the Fourier transform of the Fcrmion Green's function reads: 



I - ( I - l(Jf + + K^), H _ _L(^ + + A , : ;■: 2: ) (-)!), 



When this Green's function is analytically continued to real frequency, it is seen[146j that it has a power law 
singularity for to 2 — (uq) 2 +m 2 , and is analytic for u> below this threshold. As a result, the Fermion Green's function 
vanishes below the gap, as it would do in a superc ondu ctor |l47j and could be checked experimentally. Note however 
that the anomalous Green's function remains zero |l48j . due to the fluctuations of the phase 9 + . 
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IV. MOTT INSULATING STATE 



Till now, we have only considered the case of the continuum system Q or incommensurate filling in the lattice 
system ©. We now turn t o a lattic e system at commen surate filling. We first establish a generalization of the 
Lieb-Schultz-Mattis theorem |l49t ll5Ct Il5ll Il52l Il53l Il54j in the case of the boson-fermion mixture described by 
the Hamiltonian (J2J) . This will give us a condition for the existence of a Mott insulating state without spontaneous 
breakdown of translational invariance. Then, we will discuss using bosonization t he pr operties of the Mott state. We 
note that Mott states have been studied in the boson-fermion model in Refs. [7^. Il04j | , but not in a one-dimensional 
case. Finally, we will consider the case when the molecules or the atoms can form a Mott insulating case in the 
absence of boson-fermion conversion, and we will show that this Mott state is unstable. 



A. Generalized Lieb-Schultz-Mattis theorem 



A generalized Lieb-Schultz-Matt is theorem can be proven for the boson fermion mixture described by the lattice 
Hamiltonian Let us introduce the operator: 



U 



exp 



N 



N 



such that U^HbfU = -ffb/. Following the arguments in Ref. [I53j . one has 



(0\rfHU-H\0) 



(100) 



(101) 
(102) 



where |0) is the ground state of the system, T is the translation operator and H is the full Hamiltonian. The quantity 
v is defined by: 



1 N 1 



Nf) 



(103) 



For noninteger v, it results from the analysis of 153] that there is a state U\0) of momentum 2m/ ^ 0[2tt] which is 
orthogonal to the ground state |0) and is only 0(1/N) above the ground state. This implies either a ground state 
degeneracy (associated with a spontaneous breaking of translational symmetry) or the existence of gapless excitations 
(if the spontaneous translational symmetry is unbroken and the ground state is unique). For integer i/, the ground 
state and the state U\0) have the same momentum. In that case, a gapped state without degeneracy can be obtained. 
This state is analogous to the Mott insulating state in the half-filled Hubbard model in one- dimension ( 1 5 5| or the Mott 
insulating state in the Bose-Hubbard model with one boson per site . We note that for A = in the Hamiltonian 
PJl fermions and bosons are separately conserved, and the corresponding Fermi wavevectors are: 



k B 
kp 



N 
irNf 

~2N 



(104) 
(105) 



The momentum of the state £/|0) is thus equal to 4(fcs + kp). 
Hubbard model, 4kp — 2n is thus generalized to 4(fc^ + kp) = 



The condition to have a Mott insulating state in the 
2tt i.e. 2N b + N f = N. 



B. umklapp term 

In this Section, we provide a derivation of the umklapp term valid in the case of the lattice system (|2J). Let us 
consider the 2kp and 2ks components of the atom and molecule charge density, given respectively by Eq. (|3H|) and 
Eq. H31(l . These terms yield an interaction of the form: 

C j dxcos(2(/)b — 2ksx) cos(v / 2</> P — 2kpx) cos 
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C f 

— j cos[2cf>b + V2(j) p — 2(ks + kp)x] cos V2(j)a- 

— / cos[20t — V2<j)p — 2(fcs — kp)x] cos V2(j) a (106) 



In Eq. H106fl . the last line is the backscattering term of 134fl . and the second line is the umklapp term. Let us 
consider a case with kp ^ ks, and let us concentrate on the effect of the umklapp term. Using the rotation l|37|) . we 
can reexpress it as: 

C f 

= — j cos[v / 6(/) + — 2(ks + kp)x] cos \[2§ a 

(107) 

In the following we will consider the cases corresponding to one or two atoms per site. 



1. Mott insulating state with one atom per site 

Let us consider first the case of (ks + kp) = Then, the term (|107f) is oscillating. In second order perturbation 
theory, it gives rise to the umklapp term: 

H u^k. = j^Ti J ^cosV240+. (108) 

The condition for the appearance of the umklapp term (|108|) can be seen to correspond to having one fermion atom 
per site of the atomic lattice. Let us briefly mention two alternative derivations of (|108|l . A simple derivation can be 
obtained by considering the combination of the Aks term in the boson density with the 4fcp term in the fermion density 
in Haldane's expansion |52j. A second der ivation ca n be obtained by considering the effect of a translation by one 
lattice parameter on the phases 4> p and </>& |l53t Il56| . The expressions of the densities (|3U[I and 131|) imply that upon 
a translation by a single site <j> p — * <j) p — \f2kpa. and 4>b — > (f>b — ksa. Therefore, the combination \/6(f>+ = 2(f>b + \f2(j) p 
transforms as :v / 6</>+ — > V&(f>+ — 2{ks + kp)a. For 2(kp + ks) = ft/a, the term cos2\/60+ is invariant upon translation, 
thus leading again to flU^I . The presence of the umklapp term IjlOSJI in the Hamiltonian can result in the opening 
of a charge gap and the formation of a Mott insulating state. Since the umklapp term is of dimension 6K + this 
implies that a Mott insulating state is possible only for K + < 1/3 i.e. very strong repulsion. For free fermions, 
the Mott transition would occur at K p = 1 i.e. for weakly repulsive interaction. Thus, we see that the Josephson 
coupling is very effective in destabilizing the Mott state. In the Mott insulating state, the superfluid fluctuations 
become short ranged. Since CDW fluctuations are also suppressed, the system shows some analogy with the Haldane 
gapped phase of spin-1 chains [l57| in that it is totally quantum disordered . In fact, this analogy can be strengthened 
by exhibiting an analog of the VBS (valence bond solid) order parameter [153, [l59|. I n Haldane gapped chains, this 
nonlocal order parameter measures a hidden long range order in the system associated with the breakdown of a hidden 
discrete symmetry in the system. The equivalent nonlocal order parameter for the atom-molecule system is discussed 
in Appendix IbI 



2. Mott insulating state with two atoms per site 



Another commensurate filling, where a Mott insulating state is possible is obtained for (kp + ks) — n/a. This case 
corresponds to having one molecule (or two atoms) per site of the optical lattice. In that case, the term in 11071) is 
non-oscillating, and it gives rise to an umklapp term of the form: 

Hllk = 7^)2 / dx cos cos V2> CT . (109) 

We notice that this umklapp term is compatible with the spin gap induced by the Josephson term l|35() . When the 
Josephson coupling is large, we can make cosv^cr — * (cos ^/2~4><j) and we see that the term l|109[) becomes relevant 
for K + = 4/3. For weaker Josephson coupling, the dimension becomes 1/2 + 3/2K + , and this term is relevant only 
for K + < 1. Since K + = 1 corresponds to hard core bosons, this means that for weak Josephson coupling, the Mott 
state with a single boson per site becomes trivial. Interestingly, we note that increasing the Josephson coupling is 
enhancing the tendency of the system to enter a Mott insulating state as a result of the formation of a spin gap. If 
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we compare with a system of bosons at commensurate filling, we note however that the Mott transition would obtain 
for Kb = 2 [53. Therefore, the Josephson coupling still appears to weaken the tendency to form a Mott insulating 
state. Such tendency was also observed in [Tflj . 



C. Commensurate filling of the atomic or molecular subsystem 

When the atomic subsystem is at commensurate filling (4fcp = ^L), an umklapp term: 

^^cosVS^, (110) 

must be added to the Hamiltonian. Such umklapp term can create a gap in the density excitations of the unpaired 
atoms. However, we must also take into account the term l|35f) . This term is ordering 0- and thus competes 
with the umklapp term i|110|) . To understand what happens when 9- is locked, it is convenient to rewrite the 
umklapp term l|ll()|l as cx cos -\/8/3("\/2<^+ — </>_). The terms generated by the renormalization group are of the form 
cos ny / 8/3(v / 20+ — 0-), with n an integer. When 0_ is locked, replacing the terms e 1 ^ - by their expectation values, 
we find that all these terms vanish. Therefore, no term cos (3(ft+ can appear in the low energy Hamiltonian. A more 
formal justification of the absence of the cos (3<ft+ term in the low energy Hamiltonian can be given by noting that 
when the Hamiltonian is expressed in terms of <f>± it has a continuous symmetry <ft + — > <ft + + a and 0_ — ► 0_ + y/2a. 
As a result, terms of the form cos {3<p+ are forbidden by such symmetry. The consequence of the absence of cos/3<^+ 
terms in the Hamiltonian when 9_ is locked is that, even if the unpaired atom density is at a commensurate filling, 
the umklapp terms do not destabilize the coupled phase. However, in the opposite case of a strong umklapp term 
and a weak boson-fermion conversion term, it is the field <ft p that will be ordered. The previous arguments can be 
reversed and show that the formation of a Mott gap for the fermions will prevent the formation of the coupled phase. 
Using the method of Ref. |160l | , one can show that the phase transition between the coupled and the decoupled state is 
identical to the phase transition that occurs in two non-equivalent coupled two-dimensional XY models. This phase 
transition was studied by the renormalization group in |l6ll Il62j| . It was found that in the case of interest to us, this 
phase transition was in the Ising universality class. Thus, one expects a quantum Ising phase transition between the 
state where the fermions are decoupled from the bosons and form a Mott insulator and the state where the fermions 
and bosons are coupled and form a superfluid. 

Of course, the same arguments can also be applied to the bosons at commensurate filling, the role of the fields <ftb 
and <fip being simply reversed. 



V. RELATION WITH EXPERIMENTS 

A. Without a potential along the tubes 

To connect experiments in quasi-one-dimensional confining waveguides with theoretical models in ID, it is necessary 
to obtain estimates of the parameters that enter the Hamiltonians J2Jj0J an d the bosonized Hamiltonian H25|) . l|35|l and 
(|32|l . Since the parameters in the Hamiltonian (0) depend on the periodic optical trapping potential, we will mainly 
focus on the parameters that enter in the continuum Hamiltonian Eq. Q and in the bosonized Hamiltonian, i.e. the 
Luttinger exponent K pi the velocity, and the fermion-boson coupling A. Before giving an estimate of the parameters 
we need first to remind that, at the two-body level, there is a connection between the ID boson-fermion model 
and the quasi-lD single channel model [74|. thus we will use one or the other depending on the physical parameter 
wherein we are interested. Experimentally, molecules have been formed from fermionic atoms 6 Li|lll Il3l 111 . [2(j and 
40 K^, 0, [6(|. For 6 Li, the mass is roughly 6 times the mass of the proton, mi?( 6 Li) = 9.6 x 10 _27 kg, and for 
40 K it is mF( 40 K) = 6.4 x 10 _26 kg. Then, we need to determine the effective interaction of atoms under cylindrical 
confinement. For the interaction, we will assume a contact form, i.e. V(x — x') — gioS(x — x') and the Hamiltonian 
describing the confined atoms reads: 

H = J dr^(r,f)AVvM) + y / dr^(r,i)(wiri+^z 2 )^(r,i) 

+ i J drdrVi(r,t)Vi(r',i)f/(r-r')^(r',^(r,i), (HI) 
where r — (r j_ , z) , the second term represents the harmonic confinement potential, and: 



U(r) = g 3D 6(r), 



(112) 
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is th e ato m-atom repulsion. The coupling constant g^o is expressed as a function of the atom-atom scattering length 
a s as |l47j 

4:nh 2 a s 

93D = ■ (113) 

m 

We introduce the following decomposition of the fermion annihilation operator: 

V>(r, t)=J2 0m(r±)^ m a(z, t), (114) 

m 

where {ipmaiz), VWcr' (%')} = S mm 'S r7 . a >S(z — z 1 ), and the eigenstates of the transverse Hamiltonian <f> m satisfy the 
Schrodinger equation: 

ft 2 , mwi ,\ , , > ,1s,,, , 

~ 2^ A± + J Mr±) = ° ( " + 2 )0n(r±) ' (115) 

and the orthonormal condition J dr±(f> n (r±)(j) m (r±) = S nim . We will assume that we have a very elongated trap, with 
lu z <C uj± so that we can neglect the longitudinal confinement. The first line of the Hamiltonian (|lllf> can then be 
rewritten as: 



/h 2 
dz -—^ a {z,t)A^ a (z 7 t) + hLu (n + l/2)^l a (z 7 t)^ a {z,t) 



(116) 



If the transverse zero-point energy is much higher than the interaction energy per atom, the transverse motion is 
frozen in the ground state. It is known that virtual transitions to higher states can give rise to a divergence in ID 
sca tteri ng length known as confinement induced resonance (CIR) which is a kind of Fano-Feshbach resonance |b5 l64l 
IfiH Il63| . We will ignore CIR for the moment and restrict ourselves to consider the lowest energy level n = for the 
transverse motion. The interaction in (|lllfl can be rewritten as: 

H int = \\ 'dz J dz'^l a (z)^l a (z')V{z - z')VV Wo, CT W, (117) 

where the effective potential V reads: 

V(z-z') = y"dr ± dr' ± |0 o (r ± )| 2 |0o(r'±)| 2 C/(r ± -r' ± ,z-z'), (118) 

Using the expression of the ground state wave function of the transverse motion: 



Mrx) = ^e~^, (119) 

substituting it into l|118|) . using the definition of the interaction l|112fl . and integ rating over the transverse coordinate, 
we obtain V(z) — gi£>5(z) with the effective one-dimensional coupling |l63L Il64| : 

g 3D muj ± 

giD = — ^— = 2ha s ujj_. (120 

Eq. (|120f) can be generalized to the case in which atoms of different species, with different trapping frequencies uj±.i 
and loj_^2 are interacting with each other. In that case, 

gi D = Aha 12 -, ; r, (121) 

[U)_L1 + Ld±2) 

where an is the atom-atom scattering length. Knowing giD, we can obtain the Luttinger exponent of fermions in 
(123) as: 

<122> 

In the case of bosons, the Luttinger exponent Kb must be extracted from the Lieb-Liniger equations |53l l88l l03| . 
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Having obtained the form of the effective inte ract ions, we turn to the determination of the Josephson coupling A in 
the boson-fermion model Q). In the 3D case[6 Hl65| . the boson- fermion conversion factor is given by: 

X 3D J d 3 r^(r)V^ T (r)Vl(r), (123) 

with: 



A3, = J 4 ^ , (124) 
V m 

where a bg is the atom-atom scattering length far from resonance, AB is the width of the resonance and A/i is the 
difference of magnetic moment between atom and molecule. Using the projection on the lowest level, we obtain: 

A 2 = 2huj ± a bg AfiAB. (125) 

Knowing the interaction A in terms of the relevant physical parameters, we finally have to determine the spatial cutoff 
a to use in the bosonization and comment on the validity of bosonization approach. In the case of the optical lattice 
Eq. @, the obvious spatial cutoff is the lattice spacing. The cutoff to use in bosonization for the continuum case of 
Eq. Q is obtained in the following way. Bosonization is applicable as long as the kinetic energy of longitudinal motion 
of the particles is much smaller than the trapping energy hujj_. Thus, we have to impose the condition: HvfA ~ 
where A is the momentum cutoff 164]. The real space cutoff in the continuum case is thus a ~ A -1 ~ vf/lu±. 

The condition for perturbation theory to be valid is that the energy associated with the formation of molecules, 
XoT 1 / 2 is small with respect to the energy cutoff huj±. Therefore, perturbation theory is applicable when: 

« 1, (126) 



i.e. 

a bg (J,AB 



hvf 



< 1 (127) 



Using the values given in Refs. |l3l Il66| . we find that this parameter is small for vp 3> 3.2 x 10~ 2 m/s. Since vp 
can be expected to be of the order of lCU 3 m/s, this is not unreasonable. In fact, using the values of the trapping 
frequency given by Moritz et al.|66l| we find that: 



vp = J - (128) 



2 x 100m- 3 x 10- 34 J ■ s x 10 3 Hz 

~ V 6 x 1.6 x 10- 27 K 5 ( } 

= 4.6 x lCT 2 m/s (130) 

Therefore, we see that with 6 Li at the narrow resonance, the ratio is of order 0.7 and we can expect that our theory 
is valid qualitatively. 

Concerning K p , we find that: 

1 uj ± a bg 2ir x (69 x lO 3 ^ 1 x (80 x 0.5 x lO" 10 )™ . o . 

Kp ~ 1 -—^ 1 4.6 x 10- 2 m..s-i ' (131) 

~ 0.995 (132) 



i.e. interactions between fermions can be neglected. Since the interaction between the molecules 167J has a scattering 
length qbb — 0-6app one sees that molecules are only weakly interacting. 



B. With a potential along the tubes 



As we have seen in Sec. IV Al in the case of a two-dimensional optical lattice without periodic potential along the 
tubes, the repulsion between the bosonic molecules is weak, making the decoupling transition or the Mott transition 
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impossible to observe. To increase the effect of the repulsion, one needs to increase the effective mass of the atoms 
by adding a periodic potential along the tubes. A periodic potential can be imposed along the tubes by placing the 
atoms in a three dimensional optical lattice. The atoms experience a potential: 

V(x, y, z) = V x sin 2 (^j + V y sin 2 (^j + V z sin 2 (^j , (133) 

where A; is the wavelength of the laser radiation, and V x <C V y , V z so that the system remains quasi one-dimensional. 

The strength of the potential is measured in unit of the recoil energy — ^ {r^j as ^ = s Er,- Typical values 

for s are in the range 5 to 25. For lithium atoms |168j . the typical value of Er is 76kHz. If the potential is sufficiently 
strong, the atoms tend to localize in the lowest trap states near the minima of this potential. In our case, since the 
periodic potential along the tubes has shallower minima than in the transverse directions, the small overlap between 
the trap states in the longitudinal direction yields the single band Hamiltonian H3 ■ An expression of the 
parameters of the lattice model J2J in terms of the microscopic parameters has been derived inj^]. On the lattice, the 
Fermi velocity of the atoms and the pseudo- Fermi velocity of the molecules can be reduced by increasing the depth of 
the periodic potential in the longitudinal direction. This allows in principle to move the system near the decoupling 
transition [76j or the Mott tran sition, by reducing K p . 

A second possible setup |l6{^ | is to use a cigar shaped potential: 

V cigar (a;, y, z) = -rauj^x 1 + fi 2 r\), (134) 

with /i» 1, so that the atoms and the molecules are strongly confined in the transverse direction, and to apply a 
periodic potential: 

Aperiodic = Vq sin 2 (^pj , (135) 

in order to form the one dimensional structure described by the model (J2J). 

The main difficulty of experiments in optical lattices is that the reduction of the bandwidth results in a reduction 
of the Fermi velocity vf- Since the perturbative regime is defined by Aa 1 / 2 <C vf, this implies that by increasing 
the depth of the potential in the longitudinal direction one is also pushing the system in the regime where the 
boson- fermion conversion term must be treated non-r)erturbativelv|74j . However, in that regime there isn't anymore 
coexistence of atoms and molecules and the decoupling transition does not exist. Moreover, in that regime, the Mott 
transition becomes the usual purely fermionic or purely bosonic Mott transition. |53| 



VI. CONCLUSIONS 



We have studied a one-dimensional version of the boson-fermion model using the bosonization technique. We have 
found that at low-energy the system is described by two Josephson coupled Luttinger liquids corresponding to the 
paired atomic and molecular superfluids. Due to the relevance of the Josephson coupling for not too strong repulsion, 
the order parameters for the Bose condensation and fermion superfluidity become identical, while a spin gap and 
a gap against the formation of phase slips are formed. As a result of these gaps, we have found that the charge 
density wave correlations decay exponentially, differently from the phases where only bosons or only fermions are 
present [3- We have discussed the application of a magnetic field that results in a loss of coherence between 
the bosons and the fermion and the disappearance of the gap, while changing the detuning has no effect on the 
existence of the gaps until either the fermion or the boson density is reduced to zero. We have discussed the effect 
of a backscattering term which induces mutual locking of the density of bosons and fermions favoring charge density 
wave fluctuations resulting in a quantum Ising phase transition between the density wave phase and the superfluid 
phase. We have found a Luther-Emery point where the phase slips and the spin excitations can be described in terms 
of pseudofermions. For this special point in the parameter space, we have derived closed form expressions of the 
density-density correlations and the spectral functions. The spectral functions of the fermions are gapped, whereas 
the spectral functions of the bosons remain gapless but with an enhanced divergence for momentum close to zero. 
Finally, we have discussed the formation of a Mott insulating state in a periodic potential at commensurate filling. 
We have first established a generalization of the Lieb-Schulz-Mattis theorem, giving the condition for the existence 
of a Mott-insulating state without spontaneous breakdown of translational invariance. Then, we have discussed the 
properties of the Mott-state in the case of one atom or two atoms per site showing that in the first case the Josephson 
coupling is very effective in destabilizing the Mott state. Finally, we have considered the case when the atoms or the 
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molecules can form a Mott state in absence of boson-fermion conversion and shown that this Mott-state is unstable. 
To connect our results with experiments in quasi-one-dimcnsional confining waveguides we have derived estimates of 
the parameters that enter the bosonized Hamiltonian, as the Luttinger exponents, using the values of the trapping 
frequency and density used in experiments. We have seen that bosons are only weakly interacting and the necessary 
small fermionic Luttinger parameter required to realize a strongly interacting system, render the Mott insulating 
and decoupled phases difficult to observe in experiments. A nontrivial challenge is the experimental realization of 
the coupled Luttinger liquids phase with parameters tunable through the exactly solvable point (the Luther-Emery 
point). We suggest that a Fano-Feshbach resonantly interacting atomic gas confined in a highly anisotropic (Id) trap 
and subject to a periodic optical potential is a promising candidate for an experimental measurement of the physical 
quantities (correlation functions) discussed here. Finally we would like to comment on the fact that an interesting 
edge states physics is expected when open boundary conditions (or a cut one-dimensional boson-fermion system) are 
considered. The existence of edge states at the end of the system could lead to significant contribution to the density 
profile that could be tested in experiments. The physics of the edge states will be similar to the one of Haldane gap 
systems, like the valence bond solid model, and a study along this direction is in progress. 



APPENDIX A: CALCULATION OF THE INTEGRALS IN EQS. (JUJl- l|82jl 

In this appendix we will derive a slightly more general integral than those of Eqs. 1)81(1 and H82[). Namely, we will 
consider: 

/>oo 

9(y)= / K^u)K u (u)J x (yu)u a du. (Al) 



To find (|A1|) explicitly, we use the series expansion of the Bessel function J\ from |146| [Eq. (9.1.10)]. We find that 



k=0 



(y)=(») y[-!L) -4 / KJu)KJu)u 2k+x+a du. (A2) 

yj > \9. J A) T(k + l)r(k + A + 1) J M ; v ; 



The integral that appears in the expansion in powers of y 2 is a well known Weber-Schaefheitlin integral [l43l Il44i ll4E 
with two modified Bessel functions. Its expression is: 



K M K A u )u ^ d u = r(2fc + a + A + 1) r [k + 1 ] r ( * 



xr(fc+ 1 ' i/ + M + Q + A Vffc+ 1 ' t/ ' /i + a + V ), (A3) 



The resu lting expression of g(y) can be rearranged using the duplication formula for the Gamma function, Eq. 
(6.1.18) in |l46j |. The final expression of g is: 

2 

r (k + i + aS] r (k + s^J) r(fc + a +T) *i 



(y y ^ ~ r (k + r (fc + ktggg ) x r (fc + ^jg±A ) r (fc , 2 ; , ( „ 



fc=0 



This series expansion is readily identified with the definition of the generalized hypergeometric function 4F3 given in 
133]. So we find finally that: 



9(v) = -r o 



^ p I l+u+11+a+X ^ p I l + tz-p+q+A ^ p f l-u+^+a+X } p f l-v-y+a+X 



4 \2J r (1 + s±A) r r(A + 1) 

_ /1 + a + A + ^ + u l + a + A + i/ — u l + a + A — v + u l + a + A — v — a a + A a + A+ l 

X4-P3 -, -, -, -; 1 H , ■ 1 + A; 

\ 2 2 2 2 ' 2 2 

For v = /1, the function 4F3 reduces to a simpler 3F2 function. This leads to Eqs. l(HT|) and (|52l . 
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APPENDIX B: NON LOCAL ORDER PARAMETER FOR THE MOTT STATE 



The Mott insulating state can b e charact erized by the expectation value of a non-local order parameter as the 
Haldane gap state in a spin-1 chain [l58l Il59| . The non-local order parameter is defined as follows: 

i 

0(k, I) = {b\b k [J e-^^+E, fLh-tybi) (Bl) 
j>k 

The string operator: 

O string(M ) = ^ e -m*h+Y,jh^) (B2) 

j>k 

is a product of exponentials. As a result of its definition, we see that it is counting the number of fermions located 
between the sites k and I, either unbound or forming a molecule. To derive a bosonized expression of this operator, 
we notice that exp^inb^bj) — 1 since b^bj has only integer eigenvalues and rewrite the string operator as: 

^jWh-Y.JlM (B3) 

j>k 



Using bosonization and Eq. 1B2I) . we find: 



2ir 2 
*y (PB - Pf){x ~ x) - -j={<j)-{x) - <p-{x')) 



(B4) 



OstringC^a/) = exp 

Using (|55(1 and l|B4|l we obtain the nonlocal order parameter I|B1(I as: 

0{x,x') — (e'V? W+( x ) _ ^+( x 0)-*t( 2 p.b+/'f)(x-x')\ (B5) 

In the Mott insulating state with one fermion per site, we have 4(fcp + ks) = 2it(2pB + Pf) — 2n7r where n is an 
integer. Taking x,x' — > oo, we see that the expectation value of the order parameter is non-vanishing in the Mott 
state. 

A related VBS type order parameter can be defined as: 



0'(k,l) = ( 
In bosonized form, we have: 



E/U^l rie<*^-£-*U.*> fe/U,)> (B6) 



0'(x,x') — (e i V / |(<M :r )- , M ;i: '))+ i § (2pb+pf)(x-x')\ ^-qj^ 

and again this order parameter is non-vanishing. The physical interpretation of the non-zero expectation value of 
these nonlocal order parameters is that both bosons and fermions possess a hidden charge density wave order in the 
Mott insulator. This charge density wave is hidden as a result of the fluctuation of the density of fermions and the 
density of bosons. 
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